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ABSTRACT 



Aims. Previous studies indicate that faint and bright early-type galaxies (ETGs) present different coefficients and dispersion for their 
Kormendy relation (KR). A recently published paper states that the intrinsic dispersion of the KR depends on the magnitude range 
within which the galaxies are contained, therefore, we investigate here whether the magnitude range has also an influence over the 
values of the coefficients of the KR; a (zero point) and (slope). If the values of the KR coefficients depend on the magnitude range, 
and this fact is not considered when performing comparisons of different galaxy samples, the differences which might be found may 
be misinterpreted. 

Methods. We perform an analysis of the KR coefficients for 4 samples of galaxies, which contain an approximate total of 9400 ETGs 
in a relatively ample magnitude range (< AM > ~ 6 mag). We calculate the values of the KR coefficients in two ways: i) We consider 
the faintest galaxies in each sample and we progressively increase the width of the magnitude interval by inclusion of the brighter 
galaxies (increasing magnitude intervals), and ii) we consider narrow magnitude intervals of the same width (AM =1.0 mag) over 
the whole magnitude spectrum available (narrow magnitude intervals). We also perform simulations of the distribution of galaxies in 
the log(r e ) - < /j > e plane and compare the KR defined by the simulations with that obtained from the real galaxy samples. 
Results. The main results we find are as follows: i) In both increasing and narrow magnitude intervals the KR coefficients change 
systematically as we consider brighter galaxies, ii) non-parametric tests show that the fluctuations in the values of slope of the KR are 
not products of chance variations and that there is evidence of an underlying trend, and iii) this trend suggest a maximum of the slope 
around absolute magnitude M B 18 ± 1. 

Conclusions. The analysis of the results makes us conclude that the values of the KR coefficients depend on the width of the magnitude 
range and the brightness of galaxies within the magnitude range. This dependence is due to the fact that the distribution of galaxies 
in the log(r c )- < p > e plane depends on luminosity and that this distribution is not symmetrical, that is, the geometric shape of the 
distribution of galaxies in the log(r e )- < p > e plane plays an important role in the determination of the values of the coefficients of 
theKR. 

Key words. Galaxies: Elliptical and lenticular -Galaxies: Fundamental parameters 
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1. Introduction 

It is a well known fact that the structural parameters of 
ordinary ETGs follow the Fundamental Plane (FP) relation 
Pjorgovsky & Davis 1987] IDressler et al. 19871 1. The FP rela- 



tion is usually expressed as a correlation among the logarithm 
of the effective radius (log (r e )), effective mean surface bright- 
ness (< p > e ) and logarithm of the central velocity dispersion 
(log (ctq)), and is expressed mathematically with the following 
equation: 



log (r e ) = a log (cr ) + b < p > e + c 



(1) 



The FP relation is a direct consequence of the dynamical 
equilibrium condition (virial theorem) and of the regular be- 
haviour of both the mass-luminosity ratio and of the galactic 
structure of the ETGs along all the luminosity range. Due to its 
small intrinsic dispersion (~ 0.1 dex in r e and cr and ~ 0.1 
mag in < p > e ) the FP is considered a powerful tool in mea- 
suring galactic distances and also in studies of galactic forma- 
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tion and evolution ( Kjaergaard et al. 1993 Jorgensen et al. 1996 



Jorgensen et al. 1999||Kelson et al. 19971 1. 

A physically significant projection of the FP is the correla 
tion between log (r e ) and < p > e , known as the Kormendy rela 
tion (KR): 



<P>e 



a + p \og(r e ) 



(2) 



Several studies demonstrate that the ETGs in clusters define 
the KR with an intrinsic dispersion of approximately 0.4 mag in 
< p > e ([Hamabe & Kormendy 1987| IHoessel ef aZ. 19871 



Sandage & Peremulter 1991 



Sanda gF&Lubin 2001| 



La Barbera et al. 2003b . Some authors contend that the high 
intrinsic dispersion reported is due mainly to the fact that the KR 
does not consider the third parameter (the velocity dispersion of 
the FP) (Ziegler et al. 1999), besides the value of this dispersion 
must be slightly increased by the measurement errors as well as 
the systematic errors due to the photometry calibration and also 
due to the corrections introduced for different biases (zero point 
and color transformation, K-correction and reddening). 

Recent studies show that low-density-environment ETGs 
and isolated ETGs also follow the KR with the same coefficients 
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Table 1. Name, Number of galaxies, Magnitude range (according to the papers where the samples were taken from), 
Approximate magnitude range in the B-Filter (calculated by us), redshift and type of photometric profile (TP) (S = Sersic, 
dV = de Vaucouleurs) of the different galaxy samples compiled in this work. 



Sample 


N 


Magnitude range 


Approximate magnitude range in the B -filter 


z 


TP 


7 Abell clusters (V filter) 


626 


-16.0 >M V > -22.0 


-15.1 > M B > -21.1 


0.048 


S 


Coma cluster (Gunn r filter) 


196 


-17.0 > M Gr > -24.0 


-15.9 >M B > -22.9 


0.024 


S, dV 


Hydra cluster (Gunn r filter) 


54 


-18.0 > M Gr > -22.0 


-16.9 > M B > -20.9 


0.014 


dV 


SDSS (g* filter) 


8666 


-18.0 > M g , > -24.1 


-17.5 >M B > -23.6 


< 0.3 


dV 


SDSS (r* filter) 


8666 


-18.6 > M n > -24.7 


-17.5 > M B > -23.6 


<0.3 


dV 


SDSS (i* filter) 


8666 


-19.0 > Mi, > -25.1 


-17.5 >M B > -23.6 


<0.3 


dV 


SDSS (z* filter) 


8666 


-19.3 > M z , > -25.3 


-17.5 >M B > -23.6 


<0.3 


dV 



and intrinsic dispersion as do ETGs in clusters (Reda et al. 2004 
Nigo che-Netro et al. 2001) . On the other hand, previous stud- 
ies demonstrated that ETGs (plus bulges of spiral galaxies) 
formed two distinct families on the structural parameters plane; 
i.e. the family of bright ETGs (M B < -18) that follows 
the KR and the family of dwarf ETGs (M B > -18) with 
more disperse and heterogeneous properties (their effective pa- 
rameters range over a wide interval for the same total lumi- 
nosity) (Kormendy 1985 Capaccioli et al. 1992). Subsequently 
Graham & Guzman (2003) showed that the difference between 
the bright elliptical (bright E) and the dwarf elliptical (dE) galax- 
ies is only apparent and that there is a continuous structural 
relation between both classes. They say that the different be- 
haviour presented by the bright E and the dE galaxies in the 
\og{r e ) - < fi > e plane and the variations in the value of the 
slope of the relation do not imply a different formation mech- 
anism, rather it may be interpreted as a systematic change in 
the shape of the light-profile with galactic magnitude (M). This 
result, together with the fact that the intrinsic dispersion of the 
KR depends on the magnitude range (Nigoche-Netro et al. 2007 



Nigoche-Netro 2007) prompted us to study the behaviour of the 
coefficients of the KR as functions of the magnitude range. If in- 
deed, the values of the coefficients of the KR depend on the mag- 
nitude range within which the galaxies are contained, and this 
fact is not considered when performing comparisons of galaxy 
samples such as the dependence of the KR on the environment, 
on redshift or on wavelength, the differences which might be 
found may be misinterpreted. 

In this paper, we present a compilation of 4 samples of ETGs 
which in total contain approximately 9400 galaxies and cover a 
relatively ample magnitude range (< AM > ~ 6 mag). Using 
these data we analyze the behaviour of the coefficients of the KR 
with respect to several characteristics of the magnitude range. 
We also present simulations of the distribution of the galaxies 
on the log(r e ) - < /j. > e plane. These simulations allow us to 
reproduce in a reasonable manner the results we obtain from the 
real samples of galaxies. 

This paper is organized as follows. In section 2, we present 
the different samples used in the analysis of the KR. Section 3 
describes the fitting method used in calculating the KR coeffi- 
cients, as well as the behaviour of the /3 coefficient with respect 
to the absolute magnitude range. Section 4 presents simulations 
of the distribution of the galaxies on the \og{r e ) - < fi > e plane 
and finally in section 5 the conclusions are presented. 

2. The samples 

We use a Sloan Digital Sk y Survey (SDSS) sample of 8666 
ETGs dBernardi et al. 20 03) in filters g*, r*, i* and z* (total 
absolute magnitude range -18 > M g , > -24.1 and its equiv- 
alent in other filters) as well as a sample of 626 ETGs in the 



Johnson V filter (-1 6 > M v > -22 ) from 7 Ab ell clus- 
ters (WINGS project, IFasano et al. 20021 IVarela 20041 a sam- 
ple in filter Gunn r (-17 > M Gr > -24) with 196 ETGs from 
the Coma cluster ( Jorgensen et al. 1995] |Milvang- Jensen 1997| 
|Jorgensen 1999| Aguerri et al. 2005) and a 54 ETGs sample in 
the filter Gunn r (-18 > M Gr > -22) from the Hydra cluster 
(Milvang- Jensen 1997 1. In the 7 Abell clusters sample we in- 
clude the following clusters: A147, A168, A193, A2457, A2589, 
A2593, and A2626. All the samples are redshift-homogeneous 
(the galaxies are contained within a narrow redshift interval), 
except for the SDSS sample which cover a relatively ample red- 
shift interval (0.01 < z < 0.3). This sample is magnitude- 
limited (Bern ardi et al. 2003bl >. Besides, within large volumes 
there could be evolution effects of the parameters of the galax- 
ies. So, in order to have a representative sample of the universe 
in a given volume without any evolution effects it is important to 
consider narrow redshift intervals. Bernardi et al. 2003b recom- 
mend Az = 0.04. This value comes from the sizes of the largest 
structures in the universe seen in numerical simulations of the 
cold dark matter family of models (Colber g \et al. 200"0] l. On the 
other hand, it is also well known that the SDSS photometry un- 
derestimates the luminosity of the brightest objects in crowded 
fields (Bernardi et al. 2007 Bernardi et al. 200Tbt . To probe the 
possible evolution effects and the photometric bias of the bright- 
est galaxies, we have built a subsample from the SDSS in the 
redshift interval 0.04 < z < 0.08. This subsample has 1670 
galaxies in each filter and cover a magnitude range < AM > ~ 4 
mag (-18.5 > M g - > -22.9 and its equivalent in other filters). 
This subsample will be referred to, from now on, as homoge- 
neous sample from the SDSS. 

In Table 1 we present relevant information (number of galax- 
ies, magnitude range, redshift and type of photometric profile) 
for the samples of galaxies we use in this paper. The magnitude 
range information is given in relation to the different filters used 
(from the literature) and also the approximate range for the B- 
magnitude (calculated by us). The transformation to the B filter 
was accomplished by use of the following equations: B - Gunn 
r = 1.15 ( |Milvang- Jensen 1997|l, B - V = 0.92 (IMichard 20001 > 
and B - g* = 0.5 ( |Fukugita et al. 1996) . 

All the samples consider the photometric parameters log(r e ) 
and < yu > e corrected for different biases, such as: seeing (Saglia 
et al. 1993), galactic extinction (Schlegel et al. 1998), K correc- 
tion (Jorgensen et al. 1992; Bruzual & Chariot 2003) and cos- 
mological dimming (Jorgensen et al. 1995). These parameters as 
well as their uncertainties were taken directly from the different 
papers cited above, except for the parameters for the Abell sam- 
ples which constitute a private communication from the WINGS 
project team ( IFasano et al. 2 002 Varela2004), and those for 
Coma (data from Aguerri et al. 2005 1. In this last case, Aguerri 
et al. (2005) give information for effective surface brightness 
(fi e ) instead of effective mean surface brightness (< /i > e ), there- 
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Fig. 1. Variation of the KR slope (B) for the different samples of galaxies in increasing magnitude intervals (upper magnitude 
cut-off). This graphs shows the B coefficient obtained using the BCES Bis method. Each point corresponds to the mean value of 
the total absolute magnitude of the galaxies contained in each magnitude interval (see table A.l). Diamonds represent the SDSS 
homogeneous sample. 



fore, we transform these data following the procedure described 
by Graham & Driver (2005). Additionally, to be consistent, we 
review carefully that the photometric parameters and their errors 
were obtained by the same methods for all samples. 



It is important to note that, the photometric parameters of 
the faint and bright ETGs in the Abell sample as well as the faint 
ETGs in the Coma sample were obtained using Sersic profile- 
fits, whereas for the bright ETGs from the other samples (which 
for the SDSS sample consists of the entire sample), these param- 
eters were obtained using de Vaucouleurs r 1 ^ 4 profile-fits (see 
Table 1). In the literature we find that the use of one profile or the 
other may affect significantly the estimations of the photomet- 
ric parameters (Caon et al. 1993; Fritz et al. 2005), this in turn 
will also affect the estimations of the KR parameters (Ziegler et 
al. 1999; Kelson et al. 2000), however, we also find that the de 
Vaucouleurs profile-fits represent a good approximation to the 
Sersic profile-fits for bright galaxies (Prugniel & Siemen 1997). 
In section 3.3 it will be shown that the use of the de Vaucouleurs 
profiles as approximations to the Sersic profiles for bright ETGs 
does not produce important biases to the results obtained in this 
paper. 



Finally, an important characteristic of SO galaxies is that, in 
general, the structural properties of their bulges show approxi- 
mately the same homogeneity as those of E galaxies. Due to this 
fact, and because KR uses effective parameters, we use the pho- 
tometric parameters from the whole galaxy in the case of bright 
Es, and from the bulge in the case of bright SOs. For dwarf galax- 
ies we follow the same procedure: information from the whole 
galaxy for the dEs and from the bulge for the dSOs. 



3. Kormendy relation 

3.1. Calculation of the Kormendy relation 

The estimation of the KR coefficients may be severely affected 
by the fitting method and by the choice of dependent vari- 
able. The biases may be larger if there are measurement er- 
rors in the variables, if these errors are correlated and/or if 
there is intrinsic dispersion. The Bivariate Correlated Errors 
and Intrinsic Scatter bisector (BCES Bis) fit (llsobe et al. 1990 
Akritas & Bershad y 1996] > is a statistical model that takes into 
account the different sources of bias mentioned above, which are 
precisely those that affect our samples. In this work we use the 
BCES Bis method for the determination of the KR coefficients. 

From the photometric parameters of the different galaxy 
samples, we calculate the coefficients of the KR in different mag- 
nitude ranges, the uncertainties of the coefficients, the correla- 
tion coefficient (Pearson Statistics) and the intrinsic dispersion 
of the KR (subtracting in quadrature from the intrinsic dispersion 
in < yu > e the residues dispersion due to the measurement errors 
of < /u > e and log(r e )) (La Barbera ^f al. 20031 . According to La 
Barbera et al. (2003), for the calculation of the intrinsic disper- 
sion it is necessary to have the measurement errors in < fi > e and 
log(r<,), these errors come directly from the papers from which 
we take the galaxy samples (see section 2), while the errors in 
the KR coefficients were calculated by us following Akritas & 
Bershady (1996) in lcr intervals. 

The KR coefficients were calculated both in increasing mag- 
nitude intervals, as well as, in narrow magnitude intervals. This 
allows us, among other things, to characterize the behavior of the 
KR coefficients with respect to the width of the magnitude range 
and the brightness of galaxies within the magnitude range, so 
the results obtained may be utilized as reference for other stud- 
ies in which different magnitude ranges are used. In appendices 
A and B we present the results for the values of the coefficients 
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Fig. 2. Variation of the KR slope (/?) for the different samples of galaxies. The graphs show the fl coefficient values obtained with 
the BCES Bis method. Each point represent a 1 mag interval (mean value of the total absolute magnitude of the galaxies contained 
in each magnitude interval analyzed, see table A. 3). Diamonds represent the SDSS homogeneous sample. 



of the KR for the different samples of galaxies and for different 
magnitude intervals. 

Here it is necessary to mention that both the figures in the 
main section of this paper as well as those tables in appendix A, 
contain the photometric information given in the filters in which 
the samples were originally observed, however for comparison 
purposes of the values of the KR coefficients for different sam- 
ples, we present in appendix B the relevant figures (Figs. 1, 2 
and 7) and the corresponding tables with the photometric infor- 
mation given in the B -filter reference frame. 

3.2. Behaviour of the Kormendy relation coefficients with 
respect to absolute magnitude range 

From the analysis of the data for our different samples we notice 
that the intrinsic dispersion of the KR (<Tkr), as clearly stated 
by Nigoche-Netro et al. (2007), changes appreciably each time 
we include brighter galaxies in the samples (upper magnitude 
cut-off), that is to say when we use increasing magnitude inter- 
vals (see Tables A.l and B.l). We also notice that the correlation 
coefficient (R) for each fit diminishes considerably. Apart from 
these changes, we can also see changes in the coefficients of the 
KR and that these changes are larger than the associated errors 
for most of the cases (differences in the /3 coefficient may be as 
large as 65%). The distribution of the ft coefficient may be seeing 
on Fig. 1 (see also Fig. B.l). On the other hand, if we perform 
the KR analysis considering first the brightest galaxies in each 
sample and we include progressively the fainter galaxies (lower 
magnitude cut-off) (see Tables A. 2 and B.2), the behaviour of 
the KR parameters is similar to that described above: both the 
intrinsic dispersion and the KR coefficient-values change sys- 
tematically as we increase the width of the magnitude interval. 

We also perform the analysis of the data using magnitude in- 
tervals that are narrow (see Tables A. 3 and B.3), that is, consid- 
ering galaxy samples in magnitude intervals of the same width 



and progressively brighter. For this case we can also see changes 
in the coefficients of the KR, however, the changes are less pro- 
nounced but still larger than the associated errors for most of the 
cases (differences in the /3 coefficients may be as large as 48%). 
We also find that the intrinsic dispersion of the KR is relatively 
low in all cases and that the correlation coefficient is, on the av- 
erage, superior to 0.9. The variation of the /? coefficient when 
we consider narrow magnitude intervals may be seen on Fig. 2 
(see also Fig. B.2) where a dependence between the /3 coefficient 
and absolute magnitude appears to be hinted at (it is interesting 
to note that this distribution, in the case of the Coma and Abell 

clusters, seems to reach a maximum at Mb 18 + 1). However, 

it may be possible for the (3 coefficient to be constant (J3 = 5) and 
that the differences we find are the result of statistical fluctua- 
tions. 

The question that /3 be constant and equal to 5 comes from 
the definition of absolute magnitude in terms of effective radius 
and effective mean surface brightness, that is: 

M =<n> e -5 \og(R e ) - 2.5 log (2 tt) - 5 log(D) + 5(3) 

where R e [arcsec] and D [pc] represent the effective radius 
and the distance to the object in question respectively and < /j > e 
[mag /arcsec 2 ] is the effective mean surface brightness. 

If we consider the effective radius in kiloparsecs (r e ) and a 
constant magnitude, we obtain the KR as follows: 

<H> e = a + 5 log(r e ) (4) 

Which implies that 

P = 5 (5) 

However, the observational data for ft seem to move away 
from the expected value. In order to clear this point, it is neces- 
sary to apply non-parametric tests to make certain that the fluc- 
tuations in the values of the slope are not products of chance 
variations. 
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Table 2. Hypothesis Test for the evaluation of the figis data from the 
different samples of galaxies. Tests 1, 2 and 3 correspond to the mean 
value, run and chi-square tests respectively. The percentages refer to the 
level of confidence with which we can reject the null hypothesis. 



Test 1 Test 2 Test 3 



Sample (increasing magnitude intervals) 



Abell 


94 % 


81 


% 






Coma 


99% 


89 


% 






Abell+Coma+Hydra 


99% 






95 


% 


SDSS in each filter 


99% 


96 


% 






SDSS sum of 4 filters 


99% 


96 


% 


90 


% 


Abell+Coma+Hydra+SDSS sum of 4 filters 


99% 


96 


% 


95 


% 


Sample (1 mag interval) 












Abell 


84% 


81 


% 






Coma 


99% 


89 


% 






Abell+Coma+Hydra 


99% 






95 


% 


SDSS in each filter 


99% 


96 


% 






SDSS sum of 4 filters 


99% 


96 


% 


90 


% 


Abell+Coma+Hydra+SDSS sum of 4 filters 


99% 


96 


% 


95 


% 
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Fig. 3. Distribution of the Abell, Coma, Hydra and SDSS (r* filter) galaxies on the log(r e )- < y. > e plane. Each symbol represents 
a 1 mag interval. The circles represent the (-21, -22] mag interval. 



3.3. Hypothesis Tests for the evaluation of (3 data 

There are several non-parametric methods for the evaluation of 
sets of data (IBendat & Piersol 19661 1. One of the most popular 
is the chi-square test. This test measures the discrepancy be- 
tween an observed probability density and a theoretical prob- 
ability density (i.e. a normal distribution). Another important 
test, which does not a priori assume a specific distribution, is 
that known as run test. There is one particularly interesting non- 
parametric test (mean value test), which helps in finding out 
whether the data in question are distributed at random around 



a given value or whether this distribution is not a random one. 
For full details on the hypothesis tests see appendix C. 

In Table 2 we show the results of the application of the tests 
to the different galaxy samples in increasing and narrow magni- 
tude intervals. The null hypothesis of the mean value test (test 1) 
is that /3 has a normal distribution and that its mean value is 5 (we 
consider the data to have a measurement error equal to 10%), the 
null hypothesis of the run test (test 2) is that there is not an un- 
derlying trend in the (3 data and finally the null hypothesis of the 
chi-square test (test 3) is that the j3 data are random and that they 



6 



Nigoche-Netro et ah; The Kormendy relation. Dependence on the magnitude range. 




log(f e ) 

Fig. 4. Distribution of simulated galaxies on the log(r e )- < /i > e plane. The data distribution consists in parallel lines of slope 5 
which shift to brighter magnitude. 



follow a normal distribution. The percentages given in Table 2 
refer to the confidence level with which we can reject the null 
hypothesis. 

From Table 2 we can see that, on average, the null hypothesis 
may be rejected with a level of confidence of 95 %. This implies 
that there are strong reasons to believe that the mean value of (3 
is not 5, that there is an underlying trend in the values of /3 and 
that the distribution of these values is not normal. 

From the previous results there is a question that arise: Why 
is there an underlying trend in the values of (3 when we con- 
sider increasing and narrow magnitude intervals? A possible 
answer to this question is that the change in the value of the 
slope might be due to the fact that the distribution of the galax- 
ies on the log(r e ) - < fi > e plane depends on the luminosity 
(Fig. 3, see further details i n lVarela 20041 ID' Onofrio et al. 20061 
Nigo che-Netro et al. 2001) . It could also be due to the geometri- 
cal shape of the galaxy distribution on this plane. If the distribu- 
tion of galaxies takes a rectangular shape, fitting a straight line 
to these data will produce different results from a straight-line fit 
to a galaxy distribution that takes a triangular shape or any other 
shape. In the following section we present simulations of the 
galaxy distribution on the log(r e ) - < yu > e plane that elucidate 
clearly the effects that the shape of the distribution of galaxies 
on this plane has over the values of the coefficients of the KR. 

Here, it is important to note that the behaviour of the param- 
eters of the KR in the bright regime is similar for all the samples, 
in other words, the use of the de Vaucouleurs profiles as an ap- 
proximation of the Sersic profiles (for bright galaxies in some 
samples) does not affect appreciably the behaviour of the pa- 
rameters of the KR (see Fig. B.2 and Table B.3). For example, 
if we compare the data from the brightest part (Mb < -18) of 
the Abell sample (Sersic profile) with the rest of our samples (de 
Vaucouleurs profile) we find that the average difference between 
the /3 coefficients is 9%, which amount to the size of the errors. 
On the other hand, when we compare the behaviour of the (3 
coefficient for the heterogeneous SDSS sample and the homoge- 
neous SDSS sample (Figs. 1 and 2) we note that this behaviour 
is similar for both samples except for the brightest part, however, 
it is precisely for these galaxies for which the photometry bias 
could be more pronounced. So it is not possible to say that the 
differences found are produced by evolution effects. 



4. Simulations of the galaxy distribution on the 

log(r e ) - < /u > e plane 

To investigate whether the dependence of the KR coefficients on 
the magnitude range is due to the geometric shape of the distribu- 
tion of galaxies on the \og(r e ) - < fj. > e plane (geometrical effect), 
we perform simulations for each one of our samples of galaxies 
(Fig. 4). The simulations consist in giving values to log(r c ) in a 
similar radii range as that of the sample in question, fixing the 
faintest zero-point (co) of this sample (it is equivalent to fixing 
a magnitude) and fixing a slope of 5 (expected slope consider- 
ing the definition of total absolute magnitude). Then, we take 
the same range of \og(r e ), a ao slightly brighter (we take incre- 
ments of 0.1 mag) and the same slope of 5, and so on until we 
cover the whole range of brightness of the sample in question. 
The resulting data distribution consists in parallel lines of slope 
5 which shift to brighter magnitude. Once the data are gener- 
ated, we simulate the observed depopulation effect on the upper 
region of the galaxy distribution (Fig. 5). This depopulation is 
known as the exclusion zone (IBender et al. 19921 and may be 
characterized by a straight line (Line of Avoidance or LOA) that 
has a slope approximately equal to 2.7 dD'Onofrio et al. 2 006). 
Since we do not know the way in which galaxies locations in the 
lower part of the \og(r e ) - < /j. > e plane behave, we consider the 
following 3 cases in the simulation of this region: 

Case 1. It consists in setting a fixed limit-radius for all the 
magnitude intervals, so that the galaxies are contained within a 
triangle, one of whose sides is parallel to the < ju > e axis. 

Case 2. In this case we consider that the lower region of the 
diagram is limited by a 2.7-slope straight line, that is, there is 
an identical exclusion zone in this part of the diagram, as that 
observed in the upper part. Therefore the galaxies appear to be 
contained within a parallelogram and their distribution appears 
to be symmetric with respect to an axis that contains the barycen- 
ter of the galaxy distribution and is parallel to the log(r c ) axis, 
we shall call this axis from now on the .X^^-axis. The galaxy 
distribution is symmetrical under a reflection with respect to the 
^Csrig/w-axis and a 180° rotation with respect to a perpendicular 
line to the X/, n( ,/, r axis that contains the barycenter of the galaxy 
distribution (see Fig. 6). 

Case 3. We consider the lower region as limited in bright- 
ness, so galaxies are contained within a triangle with one of its 
sides parallel to the log(r e ) axis. 
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Fig. 5. Distribution of simulated SDSS galaxies (r* filter) on the log(r e )- < fi > e plane. Each symbol represents a 1 mag interval. 
The circles represent the (-21, -22] mag interval. Case 1: upper region of the diagram, limited by a straight line of slope 2.7 (LOA) 
and radius limit (log(r e ) < 1.6). Case 2: upper and lower regions of the diagram limited by a straight line of slope 2.7. Case 3: upper 
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22 - 




Fig. 6. Distribution of simulated SDSS galaxies (r* filter) on the log(r e )- < fi > e plane (Case 2). Each symbol represents a 1 mag 
interval. The solid line (X Bnf /„-axis) is a line of slope=0 that contains the barycenter of the distribution of the galaxies. 



The results of the analysis of the KR show that when we 
consider increasing magnitude intervals (Tables A. 4 and B.4), 
there are slope changes in all the cases, and these changes are 
always larger than the errors. We also find that the geometric 
shape of the distribution of galaxies changes each time we in- 
clude brighter galaxies in the samples (Fig. 5). On the other 
hand, when we consider narrow 1 -magnitude intervals (Tables 
A. 5 and B.5), there are also slope changes (except for case 2), 
however, the changes are less pronounced but still larger than the 



errors and if the magnitude intervals are progressively narrower, 
then the changes diminish considerably getting ever closer to the 
value of 5, just as it occurs for the real samples. Finally, if the 
magnitude interval is equal to 0. 1 mag then the slope is exactly 
equal to 5 since that is the way we define the samples. 

It is important to mention that when we consider narrow 
magnitude intervals we are able to reproduce, in a reasonable 
manner, the j3 coefficient variations and the underlying trend 
found for the real samples (except for case 2) (see Figs. 7 and 



8 



Nigoche-Netro et ah: The Kormendy relation. Dependence on the magnitude range. 



j5 

aq 4 



3 - 



t>5 

^ 4 



3 - 



SDSS (r*) Case 1 
♦ ♦ * . 

$ 

"1,1,1,1,1, 


SDSS (r*) Case 2 
♦ ♦*♦=•= 

" 1 , 1 , 1 , 1 , 1 , 


SDSS (r*) Case 3 
♦ * 

* 

— • 

_ I 

" 1 ■ 1 , 1 , 1 ■ 1 ■ 


-16 -18 -20 -22 -24 -2 



-16 -18 -20 -22 -24 -26 



Fig. 7. Variation of the KR slope (J3) for one of the SDSS simulations (r* filter). Graphs show the data for the /3 coefficient obtained 
by the BCES Bis method. Each point represents a 1 mag interval (mean value of the total absolute magnitude of the galaxies contained 
in each magnitude range analyzed, see table A.5). 



B.3). We also find that this trend seems to have a maximum 

around absolute magnitude Mb 18 ± 1. Finally, we find that 

the geometric shape of the distribution of galaxies changes sys- 
tematically when we consider brighter magnitude intervals (ex- 
cept for case 2) (Fig. 5). We must remind the reader that Case 2 
corresponds to a symmetrical galaxy distribution over the log(r e ) 
- < fi > e plane. 

Moreover, we find that both the zero point (a) and the in- 
trinsic dispersion (cr KR ) of KR change systematically when we 
consider brighter galaxies. This latter result confirms the depen- 
dence of the intrinsic dispersion with the magnitude range as 
reported in Nigoche-Netro et al. 2007. 

From the aforementioned, we may infer that the KR co- 
efficients and its intrinsic dispersion depend on the width and 
brightness of the magnitude range. This dependence is caused 
by a geometrical effect due to the fact that the distribution of the 
galaxies on the log(r e ) - < ju > e plane depends on luminosity 
and also to the fact that this distribution is not symmetrical, in 
other words, the geometric shape of the distribution of galaxies 
on the log(r c ) -<//><, plane changes systematically as we con- 
sider brighter galaxies, and so the values of the KR coefficients 
change too, because the fitting of a straight line to a set of data 
does not give the same result for slope and intercept values for 



data distributed with a rectangular shape as for data distributed 
with a triangular shape or, for that matter, with another geomet- 
rical shape. In this sense, any other systematic restrictions im- 
posed on a sample of galaxies, such as brightness cuts or effec- 
tive radius cuts, will cause changes in the geometric shape of the 
distribution of galaxies. The more pronounced these changes are 
made, the more pronounced will be the changes in the values of 
the KR coefficients. 



5. Conclusions 

We have compiled 4 samples of ETGs with information for their 
photometric parameters log(r e ) and < fi > e , these samples con- 
tain a total of ~ 9400 galaxies in a relatively ample magnitude 
range (< AM > ~ 6 mag). From the values of their photometric 
parameters, we have made an analysis of the behaviour of the 
coefficients and intrinsic dispersion of the KR with respect to 
several characteristics of the magnitude range within which the 
galaxies are contained. The results from this study are presented 
as follows: 

- We find that when we include in the samples galaxies which 
get progressively brighter (increasing magnitude intervals) 
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or if we consider galaxy samples in progressively brighter 
fixed-width magnitude intervals (narrow magnitude inter- 
vals), the KR coefficients change and these changes result 
to be larger than the associated errors for most of the cases. 
We also find that the distribution of the values of the ft coeffi- 
cient in narrow magnitude intervals might have a maximum 
at Mb ~ - 1 8 + 1 . We perform non-parametric tests on the 
P coefficient data and they indicate that the variations are 
real and that there is evidence of an underlying trend, that is, 
there is evidence that the /3 coefficient changes systematically 
when we consider brighter galaxies. 

- We perform simulations of the different samples of galaxies 
under study. The results of the analysis of the variation of the 
KR coefficients, both in increasing and in narrow magnitude 
intervals show that the coefficients depend on the width and 
brightness of the magnitude range and that this dependence 
comes as a result of a geometrical effect due to the fact that: 

- The distribution of galaxies on the log(r e ) - < /u > e plane 
depends on luminosity, and 

- That the geometric shape of the distribution of the galax- 
ies on this plane is not symmetrical (see Sect. 4 for full 
details). 

- Finally, simulations confirm the fact that the intrinsic disper- 
sion of the KR depends on the magnitude range, as asserted 
in Nigoche-Netro et al. (2007). 

From the previously mentioned results, it is very important 
to establish that if the magnitude range is not taken into consid- 
eration when performing comparisons of galaxy samples such 
as the dependence of the KR on the environment, on redshift 
or on wavelength, the differences which might be found may be 
misinterpreted. 
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Table A.l. KR coefficients for the different galaxy samples in increasing intervals of magnitude (upper magnitude cut-off). MI is the 
total absolute magnitude interval within which the galaxies are distributed, N is the number of galaxies in the magnitude interval, 
ctBis is the zero point of KR, j3 Bjs is the slope of KR, cr KR is the intrinsic dispersion of KR and R is the correlation coefficient of the 
fit (Pearson Statistics). 
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Table A.2. KR coefficients for the different galaxy samples in increasing intervals of magnitude (lower magnitude cut-off). MI is the 
total absolute magnitude interval within which the galaxies are distributed, N is the number of galaxies in the magnitude interval, 
asis is the zero point of KR, j3 Bjs is the slope of KR, cr KR is the intrinsic dispersion of KR and R is the correlation coefficient of the 
fit (Pearson Statistics). 
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17.896±0.012 


2.880±0.015 
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0.775 


-19 > M g , >-24 


8635 


18.972±0.012 


2.589±0.015 
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0.413 
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Table A.3. KR coefficients for the different galaxy samples in narrow 1 mag intervals. MI is the total absolute magnitude interval 
within which the galaxies are distributed, N is the number of galaxies in the magnitude interval, asts is the zero point of KR, /3 B is is 
the slope of KR, cr KR is the intrinsic dispersion of KR and R is the correlation coefficient of the fit (Pearson Statistics). 
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-16.0 > M v >-17.0 
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20.019±0.138 


2.537±0.276 
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-17.0 > M v >-18.0 


148 


21.009±0.051 


5.276±0.193 
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0.924 


-18.0 > Mcr >-19.0 


51 


19.087+0.111 


3.783±0.256 


0.623 
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-18.0 > M v >-19.0 


110 


20.223±().()82 
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0.310 


0.918 


-19.0>M Cr >-20.0 


14 


18.957±0.081 


3.243±0.090 


0.316 
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-19.0 > M v >-20.0 


125 


19.128±0.044 


4.983±0.093 


0.287 


0.929 


-20.0 > M Gr >-21.0 


27 


18.015±0.044 
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74 
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51 
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33 


18.246±0.071 


3.830±0.072 


0.146 


0.941 


-22.0 > M Gr >-23.() 


32 


16.879±0.035 


4.512±0.156 


0.221 
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4 
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11 
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74 


18.864±0.035 


4.966±0.091 


0.220 


0.948 


-20.0 > M g , >-21.0 


2339 


18.254±0.013 


4.463±0.021 


0.249 


0.918 


-20.0 > M r , >-21.0 


937 


18.036±0.014 


4.619±0.028 


0.233 


0.937 


-21.0 > M gt >-22.0 


3831 


17.590±0.014 


4.277±0.017 


0.245 


0.912 


-21.0 > M r , >-22.() 


3043 


17.366±0.011 


4.342±0.()18 


0.244 


0.914 


-22.0 > Mg, >-23.0 


1829 


16.874±0.023 


4.279±0.023 


0.223 


0.927 


-22.0 > M r , >-23.0 


3589 


16.624±0.014 


4.320±0.017 


0.242 


0.905 


-23.0 > Mg, >-24.0 


145 


16.386±0.086 


4.118±0.072 


0.161 


0.951 


-23.0 > M r , >-24.0 


984 


16.016±0.027 


4.224±0.028 


0.199 


0.934 
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36 


15.474±0.162 
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0.151 
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SDSS (i*) 












SDSS (z*) 








-19.0 > M it >-20.0 


31 


17.458±0.059 


4.717±0.209 


0.219 


0.937 


-19.5 > M z , >-20.0 


15 


17.331±0.038 


4.482±0.104 


0.116 


0.984 


-20.0 > Mi, >-21.0 


474 


17.926±0.018 


4.713±0.039 


0.235 


0.941 


-20.0 > M z , >-21.0 


195 


17.859±0.024 


4.838±0.056 


0.229 


0.986 


-21.0 > Mi, >-22.0 
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17.213±0.012 


4.481±0.021 


0.249 


0.918 


-21.0 > M z , >-22.0 
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17.065±0.013 
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16.572±0.012 
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0.241 
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-23.0 > Mi, >-24.0 
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-23.0 > M z , >-24.0 
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15.831+0.015 
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0.223 


0.977 


-24.0 > Mi, >-25.0 


152 


15.212±0.068 


4.240±().()62 


0.160 


0.956 


-24.0 > M z , >-25.0 


381 


15.308±0.035 


4.027±0.033 


0.158 


0.973 
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Table A.4. KR coefficients for the SDSS simulation in r* filter. Increasing magnitude intervals. MI is the total absolute 
magnitude interval within which the galaxies are distributed, N is the number of galaxies in the magnitude interval, a Bis is 
the zero point of KR, f} Bis is the slope of KR, cr KR is the intrinsic dispersion of KR and R is the correlation coefficient of the 
fit (Pearson Statistics). 
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18.232±0.022 
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>-22 
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-20 > M„ 


>-23 
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0.785 


0.945 
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>-23 


420 
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3.213±0.018 


0.365 
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-20 > M„ 


>-24 
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17.977±0.025 


4.442±0.027 
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0.907 


-20 > M r , 


>-24 


560 
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3.014±0.013 


0.381 


0.971 


-20 > M r , 


>-25 
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4.305±0.031 


1.030 


0.885 


-20 > M n 


>-25 
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18.311+0.017 


2.872±0.011 


0.388 


0.976 








Case 3 






















-20 > M n 
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245 
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4.819±0.025 
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18.474±0.027 


3.053±0.028 


0.640 
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Table A.5. KR coefficients for the SDSS simulation for r* filter. Narrow 1 mag intervals. MI is the total absolute magnitude 
interval within which the galaxies are distributed, N is the number of galaxies in the magnitude interval, asi S is the zero point 
of KR, fSgis is the slope of KR, cr KR is the intrinsic dispersion of KR and R is the correlation coefficient of the fit (Pearson 
Statistics). 
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4.508±0.043 


0.253 


0.960 
















-23 > Mr, 


>-24 
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16.697±0.068 


4.028±().()65 


0.218 


0.925 
















-24 > M r , 


>-25 


35 


17.391+0.168 


3.097±0.119 


0.140 


0.891 
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Appendix B: Figures and Tables with behaviour of 
the KR coefficients with respect to absolute 
magnitude range. The information presented 
here corresponds to photometry in the B filter. 
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Fig. B.l. Equivalent to Fig. 1 with photometric information in the B-filter. Variation of the KR slope (J3) for the different samples of 
galaxies in increasing magnitude intervals (see Table B.l). Diamonds represent the SDSS homogeneous sample. 
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Fig. B.2. Equivalent to Fig. 2 with photometric information in the B-filter. Variation of the KR slope (J3) for the different samples of 
galaxies. Each point represents a 1 mag interval (see Table B.3). Diamonds represent the SDSS homogeneous sample. 
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Fig. B.3. Equivalent to Fig. 7 with photometric information in the B-filter. Variation of the KR slope (J3) for one of the SDSS 
simulations. Each point represents a 1 mag interval (see Table B.5). 



Nigoche-Netro et al. : The Kormendy relation. Dependence on the magnitude range. 



17 



Table B.l. Equivalent to Table A.l with photometric information in the B-filter. KR coefficients for the different galaxy samples in 
increasing intervals of magnitude (upper magnitude cut-off). MI is the total absolute magnitude interval within which the galaxies 
are distributed, N is the number of galaxies in the magnitude interval, a Bis is the zero point of KR, j3 Bis is the slope of KR, cr KR is the 
intrinsic dispersion of KR and R is the correlation coefficient of the fit (Pearson Statistics). 
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19.131±0.038 


4.482±0.104 


0.116 


0.984 


-17.5 > Mb >-19.5 


505 


19.533±0.027 


4.540±0.062 


0.311 


0.894 


-17.7 > Mb >-19.2 


210 


19.740±0.035 


4.742±0.085 


0.299 


0.908 


-17.5 > M B >-20.5 


2805 


19.105±0.020 


3.941 ±0.036 


0.371 


0.816 


-17.7 > Mb >-20.2 


1799 


19.112±0.022 


4.259±0.042 


0.358 


0.850 


-17.5 > Mb >-21.5 


6645 


19.163+0.015 


3.047±0.022 


0.406 


0.741 


-17.7 > Mb >-21.2 
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19.064±0.017 


3.309±0.025 
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-17.5 > Mb >-22.5 
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19.308±0.013 
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0.716 


-17.7 > Mb >-22.2 
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19.237±0.013 


2.680±0.017 


0.414 


0.738 


-17.5 > Mb >-23.5 


8664 


19.343±0.012 


2.529±0.016 


0.422 


0.723 


-17.7 > Mb >-23.2 


8658 


19.304±0.012 


2.543±0.015 


0.413 


0.748 














-17.7 > Mb >-23.5 


8665 


19.307±0.012 


2.537±0.015 


0.413 


0.750 



18 



Nigoche-Netro et at: The Kormendy relation. Dependence on the magnitude range. 



Table B.2. Equivalent to Table A. 2 with photometric information in the B-filter. KR coefficients for the different galaxy samples in 
increasing intervals of magnitude (lower magnitude cut-off). MI is the total absolute magnitude interval within which the galaxies 
are distributed, N is the number of galaxies in the magnitude interval, a Bis is the zero point of KR, j3 Bis is the slope of KR, cr KR is the 
intrinsic dispersion of KR and R is the correlation coefficient of the fit (Pearson Statistics). 
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8658 


19.304±0.012 


2.543±0.015 


0.413 


0.748 
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Table B.3. Equivalent to Table A. 3 with photometric information in the B-filter. KR coefficients for the different galaxy samples in 
narrow 1 mag intervals. MI is the total absolute magnitude interval within which the galaxies are distributed, N is the number of 
galaxies in the magnitude interval, a Bis is the zero point of KR, j3 Bis is the slope of KR, cr KR is the intrinsic dispersion of KR and R 
is the correlation coefficient of the fit (Pearson Statistics). 
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Table B.4. Equivalent to Table A.4 with photometric information in the B-filter. KR coefficients for the SDSS simulation in 
r* filter. Increasing magnitude intervals. MI is the total absolute magnitude interval within which the galaxies are distributed, 
N is the number of galaxies in the magnitude interval, a Bis is the zero point of KR, fi Bis is the slope of KR, cr KR is the intrinsic 
dispersion of KR and R is the correlation coefficient of the fit (Pearson Statistics). 
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Table B.5. Equivalent to Table A. 5 with photometric information in the B-filter. KR coefficients for the SDSS simulation for 
r* filter. Narrow 1 mag intervals. MI is the total absolute magnitude interval within which the galaxies are distributed, N is 
the number of galaxies in the magnitude interval, aBis is the zero point of KR, /3bis is the slope of KR, ctkr is the intrinsic 
dispersion of KR and R is the correlation coefficient of the fit (Pearson Statistics). 
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19.687±0.015 


4.971±0.011 
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Appendix C: Hypothesis Tests. 

C.1. Mean value test 

Consider the case where some estimator is computed from a 
sample of N independent observations of a random variable x(k). 
Assume that there is reason to believe that the true parameter 
being estimated has a specific value Oq. Now, even if 0=0o, 
the sample value will probably not come out exactly equal to 
0o because of the sampling variability associated with 0. Hence 
the following questions arises. If it is hypothesized that 0=0o, 
how much difference between and 0o must occur before the 
hypothesis should be rejected as being invalid? This question can 
be answered in statistical terms by considering the probability of 
any noted difference between and 0o based upon the sampling 
distribution for If the probability of a given difference is not 
small, the difference would be accepted as a normal statistical 
variability and the hypothesis that 0=0o would be accepted. 

To clarify the general technique, assume that has a proba- 
bility density function of p(<&). Now, if a hypothesis that 0=0o is 
true, then p(&) would have a mean value of Oo. The probability 
that O would fall below the lower value $1-072 is: 



/4>i-a/: 
CO 



a 



p(<D)d<D = - 



(CI) 



The probability that would fall above the upper value O a /2 



is: 



Prob[<& > <D ff/2 ] 



r 



P (<D)ja> = - 



(C2) 



Hence the probability that would be outside the range be- 
tween 1-0-/2 and <t>Q,/2 is a. Now let a be small so that it is very 
unlikely that would fall outside the range between 0i-«/ 2 and 
tt / 2 . If a sample were collected and a value of were computed 
which in fact fell outside the range between 0i_ a / 2 and 0»/ 2 , 
there would be strong reason to question the original hypothesis 
that 0=$o since such a value for would be very unlikely if the 
hypothesis were true. Hence the hypothesis that 0=0o would 
be rejected. On the other hand, if the value for fell within the 
range 0i_q./ 2 and a / 2 , there would be not strong reason to ques- 
tion the original hypothesis. Hence the hypothesis that 0=0o 
would be accepted. 

The small probability a used for the hypothesis test is called 
the level of significance for the test (1- a is called the level of 
confidence). The range of values of for which the hypothe- 
sis will be rejected is called the region of rejection or critical 
region. The range of values of for which the hypothesis will 
be accepted is called the region of acceptance. The simple hy- 
pothesis test outlined above is called a two-side test because, if 
the hypothesis is not true, the value of could be either greater 
or less than 0o. Hence it is necessary to test for significant dif- 
ferences between and 0o in both directions. In other cases a 
one-sided test might be sufficient. 

C.2. Chi-Square test 

A special type of hypothesis test which is often used to test 
the equivalence of a probability density function for sampled 
data to some theoretical density function is called the chi-square 
goodness-of-fit test. The general procedure involves the use of 



a statistic with an approximate chi-square distribution as a mea- 
sure of the discrepancy between an observed probability density 
function and the theoretical density function. A hypothesis of 
equivalence is then tested by studying the sampling distribution 
of this statistic. 

Consider a sample of N independent observations from a ran- 
dom variable x(k) with a probability density function of p(x). Let 
the observations be grouped into K intervals, called class inter- 
vals, which together form a frequency histogram. The number of 
observations falling within the ;'th class intervals is called the ob- 
served frequency in the /th class, and will be denoted by fi. The 
number of observations which would be expected to fall within 
;th class interval if the true probability density function for x(k) 
were po(x) is called the expected frequency in the /th class in- 
terval, and will be denoted by F ( . Now, the discrepancy between 
the observed frequency and the expected frequency within each 
class interval is given by fi - Fi. To measure the total discrepancy 
for all class intervals, the squares of the discrepancies in each 
interval are summed to obtain the statistics 
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The distribution for X 2 is approximately the chi-square (x 2 ) 
distribution. The number of degrees of freedom, n, is equal a 
K minus the number of different independent linear restrictions 
imposed to the observations. For the case where the chi-squared 
goodness-of-fit test is used as a test for normality, the number of 
degrees of freedom for X 2 is n = K - 3. 

Let it be hypothesized that the variable x(k) has a probabil- 
ity density function p(x) = po(x). After grouping the sampled 
observations into K class intervals and computing the expected 
frequency for each interval assuming p(x) = po(x), compute X 2 
as indicated in Eq. B.3. Since any deviation of p(x) from po(x) 
will cause X 2 increase, a one-sided test is used. The region of 
acceptance is: 
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If the sample value of X 2 is greater than Xn-a< tne hypothesis 
that p(x) = po(x) is rejected at the a level of significance. If X 2 is 
less than or equal to x% a , the hypothesis is accepted. 

C.3. Run test 

The run-test is a statistical non-parametric test which does not 
presuppose a specific distribution for the data which will be eval- 
uated. A run is defined as a sequence of identical observations 
that are followed or preceded by a different observation or no ob- 
servation at all. The number of runs which occur in a sequence of 
observations gives an indication as to whether or not results are 
independent random observations of the same random variable. 

A run-test is used when there is reason to believe that the data 
present an underlying trend; that is, the probability of the same 
event changes from one observation to the next. The mean value 
and variance of the number of runs in a sequence that contains 
N observations is given by: 
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Nigoche-Netro et ah: The Kormendy relation. Dependence on the magnitude range. 



The hypothesis can be tested at any desired level of sig- 
nificance a by comparing the number of observed runs (r) 
to the number of runs that are contained within the interval 
(r n ;i-a/2, fn-a/i) where n - j. This number is calculated from the 
equation that represents the distribution of the number of runs in 
a sequence of N observations; therefore, if we suppose that there 
is no underlying trend in the data we analyse, and we test this hy- 
pothesis with a level of significance a, then, r in our data must 
satisfy the following condition r m \- a /2 < r < r n;a /2', if so we 
say that the data present no underlying trend with a significance 
level of 1 - a or that the significance level of an underlying trend 
is equal to a. 



